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Abstract 

We give improved upper bounds on the radius of the largest ball of separable states of an m-qubit 
system around the maximally mixed state. The ratio between the upper bound and the best known 
lower bound (Hildebrand, quant. ph/0601201) thus shrinks to a constant c = 1^34/27 w 1.122, as 
opposed to a term of order y/ m log m for the best upper bound known previously (Aubrun and 
Szarek, quant. ph/0503221). We give concrete examples of separable states on the boundary to 
entanglement which realize these upper bounds. As a by-product, we compute the radii of the 
largest balls that fit into the projective tensor product of four unit balls in R 3 and in the projective 
tensor product of an arbitrary number of unit balls in K" for n — 2,4,8. 

1 Introduction 

This work deals with balls of separable states around the maximally mixed state for a multi-partite 
quantum system consisting of qubits. The best previously known upper bound on the largest radius 
of such balls was given in [1] and is of order yjmlogm ■ 6~" 1 / 2 , where m is the number of qubits 
in the system. A lower bound proportional to 6~ m / 2 was found independently in [2] and [3]. The 
proportionality constant \/3 obtained in these papers was subsequently improved to y / 54/17 in [4]. 
In this contribution we obtain an upper bound that is proportional to 6~ m / 2 with a constant 2, 
thus eliminating the factor \/ m log m and closing the gap between the upper and lower bounds up 
to the multiplicative factor \/M/27 » 1.122 ( Corollary 13. 9|) . The upper bound in [1] was obtained 
by computing an upper bound to the volume radius of the set of separable states and is hence not 
constructive. However, for the purpose of quantum computation it is essential to be able to prepare 
entangled states, and states that are closer to the maximally mixed state are in general easier to 
prepare. It is therefore of interest to have explicit descriptions of entangled states that are close to 
the maximally mixed state. The upper bound presented in this contribution is obtained by explicitly 
constructing mixed states on the boundary to entanglement that have the distance in question to the 
maximally mixed state (Corollary 13. 8ft . For a detailed motivation we refer to [21 Section VI]. 

The main idea underlying the construction is based on using so-called real orthogonal designs known 
in coding theory [5]. The construction involves projective tensor products of unit balls in R™. These 
objects were introduced in [I], where a lower bound on the radius of the largest ball that fits into such 
a product was given. Namely, the projective tensor product of m unit balls in R™ contains a ball of 
radius at least n~v n ~ 1 '< 2 . As a by-product, we show that this lower bound is exact for n — 2, 4, 8 and 
arbitrary m (Thcorem l3.3|) . This is due to the existence in these dimensions of real orthogonal designs 
with rate 1. Further we describe those points on the boundary of the projective tensor product of 4 
unit balls in R 3 which are closest to the origin (Corollary IA.5[) . and hence obtain also the exact value 
of the radius of the largest ball in this product. 

The paper is organized as follows. In the next section we provide necessary definitions and sum- 
marize known results which are needed further on. In Section 3 we construct separable mixed states 
of m qubits that lie on the boundary to entanglement, but are closer than 2 • Q~~ m / 2 to the maximally 
mixed state. In the appendix we consider the projective tensor product of four unit balls in R 3 and 
give a description of those points that lie closest to the origin. We also provide the matrix of a mixed 
state of 4 qubits that is produced by our construction. 



2 Definitions and preliminaries 

In this section we recall the definition of separable cones and introduce some related objects. We also 
provide some simple properties of these objects. 
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For a vector space E we denote the dual vector space by E* . Let Ex,... ,E n be vector spaces. 
We denote the tensor product of these spaces by Ei ® • • • ® £? n . For the tensor product E ® • • • ® E 
of n identical spaces we write shorthand E® n . The dual space to the tensor product space (R n )® m 
is the space of m-th order tensors on R™, or otherwise spoken the space of multilinear maps T : 
R™ x • • • x W 1 — ► R. The norm ||T|| of a tensor T will denote the Frobenius norm (the square root of 
the sum of squares of the components of T in an orthonormal basis). 

2.1 Tensor products of convex bodies 

Let us introduce the following notion from [1]. 

Definition 2.1. Let B\ C Ex, . . . , B n C E n be convex bodies residing in finite-dimensional real vector 
spaces. Then their •projective tensor product B = B\ ® ■ ■ ■ ® B n C Ex ® ■ ■ ■ ® E n is defined as the convex 
hull of the set {x = xx ® • • • ® x n \ xx S Bx, . . . , x n £ B n }. 

For the projective tensor product Bx ® • • • eg) Bx of n identical bodies we write shorthand £?® n . 

Definition 2.2. Let i? C -B be a subset of a real vector space. Then the polar B° of B is defined to 
be the set {x £ E* \ \ (x, y) | < 1} for all y £ B. 

Proposition 2.3. Let B\ C E%,..., B n C E n be convex bodies and let B = Bx ® • • • ® B n 6e t/ieir 
projective tensor product. Then the polar B° is the set of all multilinear maps T : Ex X • • • X E n — > R 
smc/i t/iai |T(xi, . . . , x„)| < 1 /or oiZ X\ € B\, . . . , x n S -B„. 

Proof. Assume the notations of the proposition. By definition the scalar product between a multilinear 
map T and a product element xi®---®x n ,x\ € E\, . . . , x n € is given by the value of . . . , x n ). 
But the set of product elements {x\ <g> • • • (8 x n \ X\ G Bi 7 . . . ,x„ £ B n } is exactly the set of generators 
of the convex body B. The proof concludes with the application of the definition of the polar. □ 

Let dB denote the boundary of B. 

Lemma 2.4. Q Lemma 2.10] Let B C W 1 be a convex, compact and centrally symmetric body with 
nonempty interior and define rs = td\yi x( zqb Then (B°)° = B, B° is convex, compact and centrally 
symmetric with nonempty interior and max xe B° M = r g- Moreover, if y is a unit vector, then 
rs ■ y S dB if and only if r~f^y G B° . □ 



2.2 Tensor products of convex cones 

Definition 2.5. Let K\, . . . , K n be regular convex cones, residing in finite-dimensional real vector 
spaces Ei, ... , E n . Then an element w £ E\<g>- ■ -®E n of the tensor product space is called K\®- ■ -®K n - 
separable if it can be represented as a finite sum Ylk=i v i ® ' ' ' ® v n 01 product elements such that 
Vi £ Ki for all k = 1, . . . , N; I — 1, . . . , n. The K\ • • • ® if„-separable elements form a regular convex 
cone, the K\ ® • • • ® K n -separable cone. We denote this cone by K\ ® • • • ® 7^„. 

We will speak of the -fTi <g> • • • ® i^ n -separable cone as of the tensor product of the cones K\ , . . . , K n . 
For the tensor product K\ ® ■ ■ ■ ® K\ of n identical cones we write shorthand Kf n . 

Definition 2.6. Let E be a real vector space equipped with a scalar product (•, •) and let K C E be 
a convex cone. Then the dual cone K* is defined as the set of elements y £ E such that (x,y) > for 
all x £ K. 

Proposition 2.7. Let K\ C E%, . . . , K n C E n be convex cones and let K = K% ® • • • ® _fiT„ fee i/ie 
corresponding separable cone. Then the dual cone K* is the set of all multilinear maps T : Ex X • • • X 
E n — > R smc/i £/iaf T(xi, . . . , x„) > /or all x\ £ K±, . . . , x n £ K n . 

The proof of this proposition is analogous to the proof of Proposition ^. 31 

We now consider cones generated by convex bodies. 
Definition 2.8. Let B C R™ _1 be a convex body. Then the convex cone 

t > 0, x £ b\ c R™ 

is called the cone generated by B. 

In order to work conveniently with cones K C R" generated by convex bodies B C R™ -1 we 
introduce the following notations. The standard basis elements of R" will be denoted by eo, . . . , e n _i, 
whereas the standard basis elements of R"^ 1 will be denoted by fx, . . . , f n -x- Tensors on the space R™ 
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will have indices running from to n — 1, whereas tensors on the space M™ -1 will have indices running 
from 1 to n — 1 . 

Let B C R n_1 be a convex, compact and centrally symmetric body, and let K C R™ be the cone 
generated by B. Let T be a tensor of order m on R ra , satisfying the following conditions: 

m m 

7o...o = 1, ^...q,^ = whenever ^ a fc > 0, a fc = 0. 

fc=i fe=i 

Define the m-th order tensor T' on R n_1 by 

T > ai...a m = _ ^ai...a ra V 1 <«!,..., CK m < 71 - 1. 

Lemma 2.9. Assume above notations. Then the following assertions are equivalent. 

i) T e (K® m )*, 
a) T g (B® m )°. 

Proof, i) => ii). Suppose T G (K® m )* and let y 1 , ... , y m G B. Let , . . . , be the elements of y k , 
k = 1, . . . , m. Define vectors x k e M™, fc = 1, . . . , m elementwise by Xq = 1, x k = yf, I = 1, . . . , n — 1. 
Then x fc G X for all fc = 1, . . . , m and hence 



T(x 1 , . . . , X™) — Tq...o Xq" X™ + • ■ ■ Tu 1 ...a m X 1 OLl X™ m 

Ql- 1 Qtti — 1 

TO — 1 TO— 1 

= 1 - E •■■ E ^...^ c m >o. 

a i—l a m — 1 

But then 

to— 1 n— 1 

(T',y 1 ®...®y m )=T'(y\...,y m )= E ^.-^ C < 

CK1 = 1 Q^ — 1 

and by convexity of i?® m we get (T',y) < 1 for all y G B® m . Since _B is centrally symmetric, the 
projective tensor product _B® m is also centrally symmetric, which yields (T' , —y) < 1 for all y € B® m . 
It follows that T" is in the polar of B® m . 

ii) => i). Suppose T" G (B® m )° and let x 1 , . . . , x m € K. Let x§, ■ • • be the elements of x fc , 

fc = 1, . . . , m. If Xq = for some index k, then by definition of if we have x k = and I^x 1 , . . . , x m ) = 0. 
Now suppose that x k ^ for all fc. By the definition of K we then have Xq > for all k. Define vectors 

y 1 ,..., y m G M" _1 elementwise by = k = 1, . . . , m, 2 = 1, . . . , n — 1. Clearly these vectors are 

x 

elements of B. Therefore we have T'fy 1 , . . . , y m ) < 1 and 

TO— 1 71—1 

?V, ■ ■ ■ ,* m ) - r ...o4 -< l + E ' ' ' E ^ m 



Ql— 1 Q m = 1 

to— 1 to— 1 



i- E ••• E 



a i = l 



= ^ x o "(l-T'(y i ,...,y m ))>0. 

Thus in any case T(x 1 , . . . , x m ) > for every x 1 , . . . , x m G K. Proposition ^ . 71 completes the proof. □ 

Let V = E"r=i ' ' ' y<xi... am f ai ® • • • ® / am G (R"- 1 )^" 1 . Define an element x = e ® • • • ® 

e o + E^ii ' ' ' £a~=i y« ™. Qm e Q1 ® • • • ® e Qm G (R")®" 1 . 

Lemma 2.10. Assume above notations. Then the following assertions are equivalent. 

i) x € if® m , 

ii) y G £® m . 

Proof, i) ii). Suppose x € K® m . Then there exists an integer ./V € N and vectors x kl = 
(xq 1 , . . . , xJf_i) T G K, k = 1, . . . , m, I = 1,...,N such that x can be represented as a sum x = 

E^Li xU ® ■ ' ' ® a;m/ - We can assume without generality that Xq 1 > for all k, I, otherwise x kl = 
and the corresponding term in the representation of x is zero and can be omitted. Define vectors 
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yki g jgm 1^ _ \ _ m ^ i — 1, . . . , JV elementwise by y kl = %n a = 1, • ■ • ,n — 1. By the defi- 
nition of K we have y kl £ B for all jfe,i. Define further elements y l £ (R™ -1 )®™, J = 1,...,JV by 
y< = y 1 ' <g> • • • ® y w = E^ii • ■ 'E^Li ^...cjd ® ' ' ' ® Clearly y' e 5® m for all Z. Let us 

compute the components of the element y. By the definition of x we have 

N N N 

- V T 11 t™ 1 -Vt 11 r ml v 11 v ml - V r 1 ' r m V 

ya 1 ...a m — 2-~/ a i Q ™ — .Z^O x Wqi !/q,„ ^ »ai— Om' 

It follows that y = Ei^j xj' x^ l y l . Moreover, since the ... element of x equals 1, we have 

YliLi x o l x ™ 1 — 1 an d V i s actually a convex combination of the elements y l . Thus y £ B® m . 

ii) => i). Suppose y £ B® m . Since x depends affinely on j/, we can assume without restriction of gen- 
erality that y can be represented as a product element y = y 1 ®- ■ -<8>y m , where y fc = (y^, . . . , J/„_i) T € 5 
for all fc = 1, ...,m. Define vectors x k+ = (1,^* , . . . , 2/£_i) T , = (1, — 2/i , • ■ ■ , -y n -i) T <= 
fe = 1, . . . , m. Clearly x k+ £ K. Since B is centrally symmetric, we have also x k ~ £ K. For any 
m-tuple a = (a±,..., a m ) £ { — 1, +l} m define x(a) = x 1 (ai) ® • • • ® x m (cr m ), where x k {a^) is defined 
by 

>, X / X fe+ , dfc = +1. 

|_ X K , CTfe = -1. 

Then we have a;(cr) G if®" 1 for any er 6 { — 1, +l} m . It is not hard to see that x can then be written as 

1 x 



2m-l / j 



x(a) 



and is hence also an element of K® m . □ 



2.3 Mixed states of multi-qubit systems 

Let A(n) be the space of real skew-symmetric n x n matrices and let /„ be the n x n identity matrix. 
Denote the space ofnxn complex hermitian matrices by Tl(n) and the cone of positive semidefinite 
matrices in this space by H + (n). Then we can identify the tensor space (7i(n))® m with the space 
Ti.(n m ). The tensor product of elements in Tl(n) amounts to the Kronecker product of matrices. 

The density matrix of a mixed state in an m-qubit system is given by a unit trace matrix in H + (2 m ). 
If the trace constraint is not satisfied, we speak of an unnormalized density matrix. A state is said 
to be separable if its density matrix is contained in the separable cone H+(2)® m . An orthonormal 
basis of the space H(2) is given by {-^oo, ^ci, 2 a2 ' ^T a 3}> wnere °"2j 03 ar e the Pauli matrices 
and (To = ^2- By virtue of this basis we can define an isometry X : R 4 — > Ti.(2), having the values 
2^(efc) = 2 <7fc ' ^ = 0,1,2,3 on the basis vectors of R 4 and continued by linearity. The preimage 
T~ 1 [H+(2)] of the positive semidefinite matrix cone under this isometry is the 4-dimensional Lorentz 
cone Li, which in turn is generated by the unit ball BcR 3 according to Definition 12.81 



3 Mixed states of m qubits 

In this section we construct mixed states of m qubits which lie on the boundary to entanglement and 
have a distance of less than 2 • Q~ m / 2 to the maximally mixed state. 

Let C be an n-dimensional linear subspace of the space M(N, N) of real N x N matrices consisting 
of similarities, i.e. for any nonzero matrix A £ £ there exists a constant a(A) such that a(A) ■ A is 
an orthogonal matrix. Let {Ui, . . . , U n } be an orthogonal basis of C consisting of orthogonal matrices. 
Define the linear map U : R" — > L by U(v) = J\ =1 t'kUk for any vector v — (vi, . . . , v n ) T £ R™. Then 
U(v) is orthogonal if and only if ||w|| = 1. Let 0(n) be the group of orthogonal n x n matrices. The 
following concept is known from coding theory [5]. 

Definition 3.1. Let L C M(N, N) be a linear subspace consisting of similarities and let {U\, . . . , U n } C 
O(N) be an orthogonal basis of £ consisting of orthogonal matrices. If the elements of the matrices 
Ui, . . . , U n are in the set {—1, 0, +1}, then we call the set of matrices U\, . . . , U n a real orthogonal 
design. The ratio n/N is called the rate of the design. 

Let us now fix two unit length vectors Wi,Wf £ M. N and a number m £ . Let U± ,...,[/„ be as 
above and consider the multilinear map 

T m :(R") m ^R, T m : (v 1 ,... 7 v m )^wfU(v 1 ) U{v m )w } . (1) 
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Lemma 3.2. Assume above notations and let B C R™ be the unit ball. The tensor T m defined in (Qp 
is in (B® m )° . If U\, . . . ,U n form a real orthogonal design and the elements of the vectors Wi, Wf are 
in the set {-1,0, +1}, then T m /\\T m \\ 2 £ dB® m . If in addition the rate of the design equals 1, then 
\\T m \\ 2 = n™" 1 . 

Proof. Assume the notations of the lemma. Let v 1 , . . . , v m £ B be vectors. Then for every k — 1, . . . , to 
the matrix U(v k ) has singular values ||i> fc ||, . . . , \\v k \\ < 1. Therefore 



|rV,...,*nil<IKIh \\{\\A\ -IKII<i- 



vfc=l 



By Proposition 12 . 31 we then have T m £ (B® m )°, which proves the first assertion of the lemma. 

Let us prove the second one. If U\, . . . , U n define a real orthogonal design and the elements of the 
vectors Wi,Wf are in the set {— 1,0,+1}, then the components of the tensor T m can equal only —1,0 
or +1. Let r\ be the number of nonzero components. Then we have HT" 1 !! 2 = r\ and j 1 ™/]^™!! 2 has rj 
non-zero components, all of which have absolute value 1/rj. Hence the tensor J^/HT" 1 !! 2 is an element 
of the L\ unit ball in (R™)® m . But this ball is a subset of B® m , because its extreme points are pure 
tensor products of unit length vectors. On the other hand, (T m , c • T m /| |T m | | 2 ) = c for any number 
cel. Hence c-T m /\\T m \\ 2 B® m for any c > 1. Thus T m /\\T m \\ 2 must lie on the boundary of B® m . 

Finally, assume that the rate of the design equals 1, i.e. N — n. Any of the matrices Uk, k = 1, . . . , n, 
has n nonzero entries, otherwise it would not be orthogonal. Therefore out of the n 3 entries of the 
matrices Uk exactly n 2 are nonzero. Let us fix two indices ki,k 2 £ {L ■■■,n} and consider the 
corresponding basis vectors e^e^ £ R™. For any unit length vector v — (vi, . . . ,v n ) T £ R™ we have 
I l e fci (Sfe=i v kUk)ek 2 1 1 < 1> because ^fe=i v kUk is an orthogonal matrix. Hence there exists at most 
one index k 3 £ {1, . . . , n} such that the (fci, k 2 ) element of Uk 3 equals —1 or +1. But in order to have 
n 2 nonzero entries in total, there must exist exactly 1 such index. 

Let us now count the nonzero components of the tensor T m . We have 



T™... am = < U ai U am w f . 

For any (m— l)-tuple of indices (a 2 , ■ ■ ■ , a m ), the vector w = U a2 U am Wf has unit length and entries 

— 1,0 or 1. The same holds by assumption for the vector Wi. Hence these vectors are proportional to 
basis vectors of R™, say ek 2 and . By the above, there exists exactly one index k^ £ {1, . . . , n} such 
that the (fci, A^) element of Uk 3 is nonzero. Therefore T™ a = for all otx ^ k% and |T™ a \ = 1 
for ct\ = k%. Now note that there are n™^ 1 different (to — l)-tuples of indices (02, . . . ,a m ). Thus 
exactly n™ -1 components of T m are nonzero and ||T m || 2 = n m_1 . □ 

Recall that for any tensor T m £ (B® m )°, Proposition ^. 31 and Lemma [2Tl yield an upper bound on 
the radius of the largest ball contained in B® m , namely l/||T m ||. On the other hand, 

n -(m-i)/2 ig a 

lower bound on this radius [IJ Lemma 4] . From the theory of Radon-Hurwitz families it is known that 
real orthogonal designs with rate 1 do exist for the values n = 1, 2, 4, 8. We obtain the following result. 

Theorem 3.3. Let B C R™ be the unit ball, and let n have one of the values 1,2,4,8. Then for any 
to £ N +; the largest ball contained in B® m has radius n^^" 1 ^ 1 ^ 2 . □ 



Let us now study the case n — 4. A real orthogonal design with rate 1 is defined by the matrices 





( X4 -X 3 X 2 -Xi \ 

X 3 X4 —X\ —Xi 

— X 2 X\ X4 —x 3 

\ X\ X 2 X 3 Xi J 



k = 1,2,3,4. (2) 



Lemma 3.4. Put the vectors Wi, uif equal to the basis vector e± and consider the tensor T m defined by 
{IP with design (0). Then T™ a 7^ if and only if the indices 1,2,3 appear in the sequence a\, . . . , a m 
either all an even number of times, or all an odd number of times. 

Proof. We have the relations C/ 4 = 7 4 , U k U[ = -UiU k for 1 < k ^ I < 3, U 2 — —I4, for k ^ 4, and 
U kl Uk 2 = ffca f° r an y even permutation {k\, k 2 , k 3 ) £ S 3 of (1, 2, 3). This leads to the relations 

T m = t™ 1 " 1 1 < k < m- 

- 1 ai...Q i _ 1 4a H1 ...t> m - L ai...a k _ 1 a k + 1 ...a m ! — — 

T m = _T m ~ 2 1 < k < m n =/= A- 

M a 1 ...a k - 2 uau k + 1 ...a m i oi...tt l _ i o t+ i...o m i — — ~ 

T '™...a fc _ 2Q/ 3a fc + 1 ...Q m = - T '™...a fc _ 2 /3 Q a t! + 1 ... Qm i 2<fc<TO, l<a^/3<3; 
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Here a(a/3j) is the sign of the permutation (a, /3, 7). 

Consider a component T™ a of the tensor T m . Above rules allow us to reduce this component to 
components of tensors of smaller order. Namely, if the indices 1,2,3 appear in the sequence a\, . . . , a m 
either all an even number of times, or all an odd number of times, then T™ am eventually reduces 
to ±T° = ±wfwf = ±1. If the opposite is the case, then T™ Qm eventually reduces to ±T„, where 
a G {1,2,3} is that index whose parity is different from the parities of the two other indices. But 
T 1 = e4, and hence T™ Qm = in this case. □ 

Let us now return to the tensor products of balls in R 3 . 

Corollary 3.5. Let B CM 3 be the unit ball. Let T m be the m-th order tensor on R from the previous 
lemma and define an m-th order tensor T m on R 3 by T™ Qm = T™ Qm , a\ . . . a m G {1, 2, 3}. Define 
further an m-th order tensor T m on R 3 by T™ Qm = T™ + a u\ . . . a m € {1,2,3}. Then for any 
meN we have f m ,f m G (B® m )° ; f m /\\f m \\ 2 ,f m /\\f m \\ 2 G dB® m . 

Proof. The tensors T m , T m are defined by the real orthogonal design (U\, U2, U3), where the matrices 
Uk are given by @, with Wi = e^, and Wf = e^ for T m , uif = U\e± = — e\ for T m . The corollary now 
follows from Lemma [3~2l □ 



Lemma 3.6. Assume the notations of the previous corollary. Then for even m we have ||T m || 2 = 
(3 m + 3)/4, and for odd m we have \\f m \\ 2 = (3 m + l)/4. 

Proof. Let m be even. By the definition of T m and by Lemma 13.41 we have T™ Qm ^ if and only if 
1,2,3 all appear in the sequence otx, . . . , a m an even number of times. Hence we get 



\T r < 



/ \ ™/ 2 / \ m/2-ki , 

El m \ \ I m \ \ -> / m — Zki 

{ 2ki 2k 2 2k 3 I = 2^ [ 2fci ) ^ [ 2k 

fei+fc 2 +fc 3 =m/2 v 7 fci=0 v 7 k=0 v 

1 m/2 / If I ) I 



, ^ , ,2™- 2fc - 1 = - l + -[(2 + ir + (2-lH = i (3'" + 3). 

Let now m be odd. Then T™ Qm ^ if and only if 1 appears in the sequence a\, . . . , a m an odd 
number of times and 2,3 an even number of times. Hence we get 

llT m ll 2 = \ ( m 

11 11 2^ \ 2fci + 1 2k 2 2k 3 

k 1 +k 2 + k 3 ={m-l)/2 v 

(m-l)/2 (m-l)/2-fei (m-l)/2 

\p / m \ ^ /m^fci-lNl^p/ m 1 9 m-2fc 1 
I 2fei + 1 / 2^ I 2k ) 2 + L \2ii + l 

fel=0 v 7 fc=0 v 7 fei=0 v 7 

= l|i + I[(2 + ir-(2-in| = i(3" l + l). □ 

Define y = Yl a a m =o yai...a m e ai <E> • • • <8 e Qm € (R 3 )® m , where the coefficients are given by 
2/0...0 = 1, 

whenever au > 0, TT afe = 0, 



Vai. 



k=l k=l 
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y<y\...a ri 



-T™ Qm whenever > 0, m even, 



3 m + 3 



y ai ...a m = 3m + 1 ^...a m whenever > 0, to odd, (3) 



where T m ,T m are the tensors defined in Corollary 13.51 By combining Corollary 13. 5[ Lemmata 13.61 
12.101 and using the fact that the cone L4 is generated by the unit ball B C R 3 we obtain the following 
result. 

Lemma 3.7. The element y G (R 4 )® m defined in lies on the boundary of the cone Lf m . □ 

Recall that there exists an isometry 1 between M 4 and 7i(2) that takes the cone L4 to the cone of 
positive semidefinite matrices H + {2). 
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Corollary 3.8. Let m e N+. The matrix p = jk^l®" 1 ^), where y G (R 4 )®" 1 is given by 
describes a separable mixed m-qubit state on the boundary to entanglement. If m is even, then p has a 

distance of 2~ m / 2 ^/ 3m 4 +3 from the maximally mixed state 2~ m ■ J 2 m • If m is odd, then p has a distance 
of 2~'"/ 2 ^ 3tT1 4 +1 " from the maximally mixed state. 

Proof. The first statement follows from the fact that the map 2 ^/ 2 T® m takes the cone Lf m to the 
cone H+{2)® m . 

Let us show the second statement. By definition the tensor y has a distance of ||T m || _1 from e® m 
for even m and of HT™!!" 1 for odd m. The maximally mixed state is the image of ef m under the map 
_X_j®m ^ ^his ma p reduces all distances by a factor of 2 m / 2 , because T® m is an isometry. Application 
of Lemma 13.61 completes the proof. □ 



Corollary 3.9. The radius of the largest separable ball of mixed states around the maximally mixed 
state for an m-qubit system is bounded above by 



2-6 m / 2 J Yq^ i m+1 , m even, 

2 ' 6 ~ m/2 \/l5^"> m odd. 

In any case, the radius of the largest separable ball is bounded above by 2 • 6 _m / 2 . 
In [3] we derived the following lower bound on this radius: 



r, =J—.6- m / 2 

1 lower w u 



The ratio of these bounds then obeys 



< 2 J— = J— « 1.122. 

Tlower V 54 V 27 

The radius of the largest separable ball around the maximally mixed state is thus determined up to a 
constant factor. 
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A Projective tensor products of four unit balls 

In this section we study the projective tensor product of four unit balls in R 3 . 

Let B C K 3 be the unit ball. Let us explicitly write out the components of the tensor T 4 = T € 
(B® 4 )° defined in Corollary O 

Taa/30 = 1 V a, (3 = 1, 2, 3; 

Tapa/3 = -1 V 1 < a ^ < 3; 

T a 0/3a = 1 Vl<a^/3<3; 

T a 0-yS — for all other index quadruples (a, (3, 7, 5). (4) 



By Lemma 15761 we have ||T|| = V21. We shall show that the tensor T is essentially the only element 
in (B m )° which has norm y/21. 

A similar result for (B® 3 )° was proven in 0]. Namely, Lemma \2. 41 and Theorems 3.8 and 3.9 in [1] 
yield the following result. 
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Lemma A.l. Let B C M 3 be the unit ball. Then max M / e(B ®3)o ||M'|| = \/7. For any M' G (B® 3 )° 
such that \\M'\ \ = y/7 there exist orthogonal 3x3 matrices U, V, W such that M' — (U ® V ® W)(JVf), 
where the tensor M is given by 

i; 

-1, a = 2, 3; 
-1, a = 2, 3; 
-1, a = 2, 3; 

0, /or a// oi/ier index triples. (5) 



M m = 

M laa = 

M a la = 

M aa l = 

M a Bj = 



The lemma implies that the tensor M given by (0 lies on the boundary of (B® 3 )°. The next result 
describes the face of B® 3 which is dual to M. This face is generated by the set {a; (£> y (£> z \ x, y, z E 
B, M(x, y, z) = 1}. 

Lemma A. 2. Let the tensor M be defined as in fifty. Then for x,y, z € B we have M(x,y,z) = 1 if 
and only if there exist angles £i, £2, £3, V such that £1 + £2 + £3 = and 

I cos£i \ / COS £ 2 \ / COS £3 \ 

x = cos ip sin £1 , y=\ cos ip sin £2 , z = cos <p sin £3 . (6) 
y sinyjsin^i / y sin ip sin £2 / y sin cp sin £ 3 J 



Proof. Assume the conditions of the lemma. The condition M(x,y, z) = 1 can be rewritten as 

y = i- 




Here Zk, k — 1, 2, 3 are the elements of z. The eigenvalues of the symmetric 3x3 matrix in the middle 
are Ai = \J z\ + z 2 + zf, A2 = — Ai and A3 = —z\. If Z2 = 23 = 0, then the corresponding eigenvectors 
are the basis vectors /1, /2, /3- If z 2 + £3 > 0, then the normalized eigenvectors are 



— zi — A 
% /2A(z 1 + A) 

% /2A( 2 ; 1 +A) 

\ x /2\(z 1 +\) / 



, A — Ai, A2; 





g3 



, A = A 3 
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An inspection of this eigenvalue decomposition reveals that the condition M(x,y,z) — 1 implies |x| = 
\y\ = \z\ — 1 and the last two components of the three vectors x, y, z must all be proportional. Therefore 
these vectors have the form Insertion of ^ in the equation M(x,y, z) = 1 yields 

cos£i cos £ 2 cos £3 — cos£i sin£ 2 sin £3 — cos£ 2 sin£i sin £3 — cos £3 sin£i sin£ 2 = cos(£i + £2 + £3) = 1- 

This completes the proof. □ 

Lemma A.3. Let M' £ (B® 3 )° be such that \\M'\\ = \/7 and M'(x,y, z) = whenever x,y, z 6 B 
and M(x, y, z) = 1, where M is the tensor defined by fifty. Then there exist £ S K, cr^, 7Tfe G { + 1, —1}, 
A; = 1,2,3 such that 



M' U2 = M[ 21 = M' 2ll = 
M( 13 = M[ 31 = M' mi = 
M m = M lQ/3 = < lj9 



^322 



: -M' 222 
: -^333 
= Kpi 

0-3 sin C, 



cos £, 
sinC, 

0, a,j3 = 2,3; 



<ri sinC, M232 = ^2 sinC, 

M 233 = n l cos C, M 323 = ^2 COS C, M^ 32 = 7T 3 COS C, 



3 

E 

3 

E 

fc=l 



Ok 



1: 



(7) 
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Proof. Assume the conditions of the lemma. Let Xk, Uk, Zki k — 1, 2, 3 be the components of the vectors 
x,y,z. 

Let us first define x, y, z by putting £i = and £3 = —£2 in @- Then x = f\, Z\ = yi, z 2 = —y2, 
z 3 = — y 3 and \y\ — 1. By the preceding lemma any triple (x, y, z) satisfying these relations fulfills the 
conditions x,y, z € B and M(x, y, z) — 1. Hence we have 

3 33 

M'(x, y,z) = J2 Mi al y ayi M '^y a yp = o 

a=l a—1 /3=2 

for any unit length vector ygl 3 . By homogeneity the condition |y| = 1 can be dropped and above 
quadratic form in the variables yi, 2/2)2/3 must be identically zero. This leads to the conditions 

M[ aa = 0, M[ al = M' na V a; M[ 23 = -M[ 32 . 

Similar relations on the components M' al g and Af^gj, a, (3 — 1, 2, 3 are obtained by putting £ 2 — or 
£3 = in ^ and following the same lines of reasoning. Combining these results yields the existence of 
a, 6, ci, C2, C3 £ K such that 

M' ll2 = M[ 21 = M 2n = a, 
M[ 13 = M{ 31 =M' 3n = b, 
M' laa = M' ala =M' aal = 0, a =1,2, 3; 

M' l23 = -M' l32 = ci, 

M' 2l3 = -M' 3l2 = c 2 , 

M' 231 = -M' 321 = c 3 . 

Inserting this and © in the relation M'(x,y,z) — yields 

(a cos ip + b sin <p) (cos £1 cos £ 2 sin £3 + cos £1 sin £ 2 cos £ 3 + sin ^ cos £ 2 cos £3 ) 
+ sin £1 sin £ 2 sin £3 (M 222 cos 3 (p + (M 322 + M 232 + M 223 ) cos 2 <p sin (p 
+ (M^ 33 + M' 323 + M' 332 ) cos tp sin 2 tp + M^ 33 sin 3 tp) = 0. 

This relation must be valid for all ^,^1,^2,^3 such that £1 + £2 + £3 = 0. Using this last relation, we 
obtain 

sin £1 sin £2 sin £3 [a cos <p + b sin tp + M 222 cos 3 tp + (M 322 + M232 + M223 ) cos 2 tp sin tp 

+{M 233 + M' 323 + M' 332 ) cos tp sin 2 tp + M' 333 sin 3 tp)] = 0. 

Putting the independent coefficients of this trigonometric polynomial to zero, we obtain the additional 
relations 

M' 222 = -a, 
M 322 +M 232 + M 223 = -b, 

M 233 + M 323 + M 332 = 

M' 333 = -b. (8) 

Now consider the condition M' € (B® 3 )°. This condition implies that the matrix formed of the 
elements M[ aB , a, ft = 1,2,3 has singular values not exceeding 1, otherwise we would find vectors 
y,z € B such that M'(f\, y, z) > 1. The singular values of this matrix are given by ^/ o? + b 2 + c\ 1 0. 
Here the first value has double multiplicity. In a similar fashion we can consider the matrices formed 
of the elements M' aig , M' a/31 , which yields the relations 

a 2 + b 2 + c\ < 1, a 2 + b 2 + c 2 < 1, a 2 + b 2 + c 2 < 1. (9) 

A similar argument leads to the conclusion that the vector formed of the elements M' a23 , a — 1,2,3 
must have a norm not exceeding 1. Permuting the indices, we get the inequalities 

4 + M l2 223 + M ,2 323 < 1, c\ + M'\ 32 + M' 2 32 < 1, c| + M'223 + M ,2 233 < 1, 

c 2 + M'\ 22 + M'\ 32 < 1, c 2 + M' 2 232 + M ,2 233 < 1, c 2 + M'\ 22 + M'\ 23 < 1. (10) 
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It follows that 



| \M'\ | 2 = 4(a 2 + b 2 ) + 2(c 2 + c 2 + c 2 ) + M'i 22 + M'\ 32 + M' 2 223 + M ,2 233 + M'\ 23 + M ,2 332 <a 2 +b 2 + 6. 

Finally using the condition ||M'|| 2 = 7, we get a 2 + b 2 — 1 and inequalities ©, (JTOJ) must actually 
be equalities. Hence there exists ( such that a — cos(, b — sin(, c\ = a = C3 = 0, |M3 22 | = 
1^232 1 = I ^223 1 > l M 2 33 l = l M 32 3 l = l M 3 32 l: 1^322 1 2 + I ^233 1 2 = 1- Combining this with ©, we obtain 
| A^3 22 1 = |sin£|, IM233I = I cos CI and equations ([7]) are verified. □ 

Theorem A.4. Let B C R 3 be the unit ball. Then max T , e(B »4 ) o ||T'|| = V2T- For any T E (B m )° 
such that | |T"| | = there exist orthogonal 3x3 matrices U, V, W, X such that (U ® V ® W <g> X)(T') 
can be obtained from the tensor T defined by by some permutation of indices. 

Proof. The tensor T defined by (|4]) has norm ^/Z\. Hence maxx'e(B®<t)° > v2L 
Let now T E (B m )° . Define three 3rd order tensors M \ M 2 , M 3 on M 3 by 

M k aPl = T' aPlk , a,/3, 7 ,k= 1,2,3. (11) 

Then for any u = ( Ul , u 2 , u 3 ) T £ B we have ^Li u fe M fe e (B® 3 ) . In particular, M fc e (B® 3 )° for all 
fc = 1,2,3 and by Lemma[0 ||M fc || 2 < 7. Therefore ||T'|| 2 = Y?k=i \\M k \\ 2 < 21, which proves the 
first assertion of the lemma. 

Suppose now in addition that ||T"|| = v2L Then necessarily ||A/ fe || 2 = 7 for all k = 1,2, 3. We will 
bring T' to the form (j4|) by successive application of orthogonal transformations in the factor spaces 
and permutations of the indices. 

By Lemma [A. 1 1 there exist orthogonal transformations U, V, W such that (U ® V ® W)(M 1 ) = M, 
where M is given by Let us assume without restriction of generality that M 1 = M. Let 

x,y,z E B be such that M 1 (x,y, z) = 1. Then for any u — (ui,u 2 ,us) T E B we have T'(x,y, z,u) = 
Efe=l u kM k (x, y, z) = ui+u 2 M 2 (x, y, z)+u 3 M 3 (x, y, z) < 1. Thus the vector (1, M 2 (x, y, z), M 3 (x, y, z)) T 
must have a norm not exceeding 1, and cos pM 2 {x, y, z) + sin ipM 3 (x 1 y, z) = for all y. 

Hence by Lemma [A.3I for any </j E [— ir, 7r] there exist £(99) G °~k(<p), Ttk{<p) £ { + 1- — 1}) & = 1,2,3 
such that the tensor M' — cos ipM 2 + smp>M 3 satisfies conditions (0 with £ = CI 1 / 3 ); °~k = c/c (</?), 
TTfe = nk{<p), k = 1, 2, 3. In particular, we have 



T X122 T 1123 \ f COS tp \ = f COs((p) 
T 1132 T 1133 / V Sin( ^ / V sin C(y) 



V <y9. 



Hence the linear map given by above 2x2 coefficient matrix is norm preserving and consequently 
orthogonal. Let us denote this coefficient matrix by C. Applying the orthogonal map diag(l, C) to the 
last factor space, i.e. performing the transformation T' i— > (I 3 8/38/3® diag(l, C))T' , we make this 
coefficient matrix equal to the identity matrix J 2 . 

We can hence assume without restriction of generality that T{ 122 — T[ 133 = 1, T{ 123 — T{ 132 = 0, 
and £(</?) = ip. Further we have cos <pT 3222 + shnpT 3223 — o~\ [ip) sin C,{f) — °~i if) sm f f° r a ^ V- It 
follows that o~i(ip) is constant as a function of ip. Similarly we can deduce that o~k(<p), 7Tfc(v) are constant 
for all k = 1, 2, 3. In the sequel we omit the dependence on p and write just cr^, 7Tfc. 

In (llip we defined the 3rd order tensors M k by putting the last index in T' to k. We can also 
choose any of the other three indices and follow similar lines of reasoning. This leads in addition to the 
relations ax + 02 + 03 = — 1, tt\ + 7r 2 + n 3 = —1 to the conditions m + o 2 + a 3 = — 1, ax + 7T2 + 7T3 = — 1, 
01 + 7r 2 + 0-3 = —1, 7Ti + ct 2 + 7r 3 = — 1, (Ti + £t 2 + ^3 = — 1, 7Ti + 7T 2 + cr 3 = — 1. Therefore a k = n k for 
all = 1,2,3, and two of the numbers ak equal —1, while the third one equals 1. By a permutation of 
the first three indices of T' we can achieve that ax — 1, a 2 = a 3 = — 1. Finally, if we now apply to T' 
the map diag(l, —1, —1) <g> I 3 (g> ^3 (E> ^3, followed by a permutation of the first two indices, we obtain 
the tensor T defined by Q. □ 



Lemma 12.41 now yields the following result. 

Corollary A. 5. Let B C K 3 be the unit ball. The largest ball that fits into the body i3® 4 has radius 
yl/21. Let Y E dB® A be of norm yl/21. Then there exist orthogonal 3x3 matrices U, V, W, X such 
that 21(U ®V ®W ® X)(Y) can be obtained from the tensor T defined by by some permutation of 
indices. □ 
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B Example of a mixed state of 4 qubits 



In this section we explicitly construct a mixed state of 4 qubits which has a distance of ^1/336 from 
the maximally mixed state and lies on the boundary to entanglement. 

For rn = 4 the mixed state constructed in Corollary 13.81 is given by the matrix 
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Corollary B.l. The matrix p defined above describes a mixed J^-qubit state on the boundary to entan- 
glement. Its multiple ^jj-p has a distance of yS/Tl from the identity matrix Ii&. 

Proof. The first statement follows from Corollary 13.81 The second statement is an easily verifiable 
consequence of similarity relations. □ 

This is the first explicit example of a mixed state on the boundary to entanglement whose multiple 
has a distance strictly less than 1 from the identity matrix. 
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